Junctions of Spin-Incoherent Luttinger Liquids with Ferromagnets and 

Superconductors 
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We discuss the properties of a strongly interacting spin-charge separated one dimensional system 
coupled to ferromagnets and/or superconductors. Our results are valid for arbitrary temperatures 
with respect to the spin energy, but require temperature be small compared to the charge energy. 
We focus mainly on the spin-incoherent regime where temperature is large compared to the spin 
energy, but small compared to the charge energy. In the case of a ferromagnet we study spin 
pumping and the renormalized dynamics of a precessing magnetic order parameter. We find the 
interaction-dependent temperature dependence of the spin pumping can be qualitatively different in 
the spin-incoherent regime from the Luttinger liquid regime, allowing an identification of the former. 
Likewise, the temperature dependence of the renormlized magnetization dynamics can be used to 
identify spin-incoherent physics. For the case of a spin-incoherent Luttinger liquid coupled to two 
superconductors, we compute the ac and dc Josephson current for a wire geometry in the limit of 
tunnel coupled superconductors. Both the ac and dc response contain "smoking gun" signatures that 
can be used to identify spin-incoherent behavior. Experimental requirements for the observation of 
these effects are laid out. 

PACS numbers: 71.10.Pm,73.21.-b,73.23.-b 



I. INTRODUCTION 

An exotic and distinctive feature of interacting one- 
dimensional systems is the phenomenon of spin-charge 
separation where the elementary excitations of the 
fermionic system are decoupled charge and spin bosonic 
modes that propagate at different velocities*^ 3 - At low 
energies, Luttinger liquid theory^ has been quite success- 
ful at describing the properties of one-dimensional sys- 
tems such as quantum wires^ and nanotubesi^ One of 
the elegant features of the theory is the ease with which 
one can obtain exact results, thanks to the quadratic 
form of the Hamiltonian describing the low energy prop- 
erties, 



H 



ID 



dx 
2^ 



1 



-^(d^iix))* + Ki(d x <t>i(x)y 



(1) 

where i — p,a stands for the charge and spin sectors and 
9i and <^ are bosonic fields satisfying [Oi(x), <f>j(x')] = 
—i^SijSgn(x — x'), and representing charge/spin density 
fluctuations and charge/spin current density fluctuations, 
respectively, and Vi are the collective mode velocities. 
The parameters Ki describe the interactions in the 1-d 
system. For repulsive interactions K p < 1, and for SU(2) 
symmetry in the spin sector K a — 1. Generally, the the- 
ory is applied under the assumption that both the spin 
and charge modes operate at low energy. But at low den- 
sities, n, where r s = (2nas) _1 3> 1 (with clb the Bohr ra- 
dius for the material), strong interactions open up a wide 
window between the characteristic energies of the charge 
and the spin sectors by suppressing the spin exchange 
energy E a while enhancing the charge energy E p . When 



the temperature T is raised into this window of energy, 
E a <C fcfli 1 <C E p (where Boltzmann's constant), the 
spin sector consists of thermally randomized spins that 
are no longer described by HJ while the charge sector 
effectively becomes a spinless Luttinger Liquid (LL). A 
one-dimensional system in this regime is known as a spin- 
incoherent Luttinger liquid (SILL). 9 

Several properties of the SILL have already been es- 
tablished theoretically, 9 but the experimental effort has 
been slowed due to the difficulty of reaching the afore- 
mentioned window of energy. Nevertheless, there are 
experimental indications of this regime in high quality 
quantum wires grown with the cleaved-edge technique, 10 
in recent split-gate devices^ and possibly even in low- 
density carbon nanotubes; 12 Given the recent progress in 
manufacturing hybrid devices consisting of superconduc- 
tors (SC) and ferromagnets (FM), and one-dimensional 
systems such as nanotubes, we propose here that the 
properties of the SILL be experimentally investigated by 
studying junctions of superconductors and ferromagnets 
with a SILL. As we show below, these SILL hybrid sys- 
tems lead to distinctive behaviors (relative to a LL or 
non-interacting 1-d system) that will allow the SILL to 
be identified in various types of transport measurements, 
and in the dynamics of a coupled ferromagnet. 

In early work on transport in the spin-incoherent 
regime, Matveev found that a SILL adiabatically con- 
nected to non-interacting leads resulted in a universal 
reduction of the conductance of a single-channel quan- 
tum wire from 2e 2 /h to e 2 //ti 13 i 14 This effect has been 
ascribed to the deactivation of the spin channel, whereby 
the spin excitations are reflected while the charge ones 
propagate freely. In a recent related work, the authors 
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studied the problem of a SILL adiabatically connected 
to two superconductors. It was found that the critical 
Josephson current through such a SC-SILL-SC systems 
suffers the same fate: the critical value of the Josephson 
current is halved relative to the LL caseJ^ These univer- 
sal results follow from the assumption of adiabaticity of 
the contacts between the SILL and the leads/SC, which 
may not be satisfied in some experimental situations. 



In this work we study the transport properties of a 
SILL in the opposite limit where it is contacted via tunnel 
junctions to a superconductor and/or a ferromagnet. We 
will discuss several scenarios that may be implemented 
experimentally with existing technologies. The central 
aim of this work is to contribute to our fundamental un- 
derstanding of the SILL by addressing its properties in 
hybrid structures, and to add to the arsenal of experi- 
ments that can be carried out to probe it and infer its 
existence. 



Our main results are the following. For the case of 
a FM tunnel coupled to a SILL, we compute the spin 
current pumped into the SILL due to a time-dependent 
magnetization in the ferromagnet. The algebraic form of 
the expression for the pumped spin current is the same 
as that for non-interacting fermions, or a Luttinger liq- 
uid, but the coefficients appearing have a temperature 
dependence characteristic of the SILL. We argue the spin 
transport is diffusive and under this assumption compute 
the renormalization (due to the pumped spin current and 
its back-flow) of the Gilbert damping a of the magneti- 
zation motion in the FM. This also has a characteristic 
temperature dependence in the spin-incoherent regime. 
For the case of a SC tunnel coupled to a SILL, we evalu- 
ate the ac and dc Josephson current in a wire that can be 
either side-coupled (bulk-coupled) or end-coupled to the 
SC. We find the dc Josephson current suffers an exponen- 
tial decay in real space due to the incoherence of the spin 
sector (in contrast to the inverse length suppression that 
was found for adiabatically contacted SCs 15 ). Finally, 
we show that the ac Josephson effect also contains im- 
portant information about spin-incoherent effects in its 
voltage and temperature dependence. 



This paper is organized as follows. In Sec. [TI] we review 
bosonization for a strongly interacting electron system 
for energy scales small compared to the charge energy, 
but arbitrary compare to the spin energy. In Sec. IIIII 
we study spin pumping and magnetization dynamics in 
the FM-SILL system. In Sec. |IV] wc study ac and dc 
Josephson effects in the SC-SILL-SC system. Finally in 
Sec. [V]we discuss prospects for experimental realizations 
of the physics we discussed and directions for future work. 
Some technical details and results appear in various ap- 
pendicies. 



II. REVIEW OF BOSONIZATION FOR 
STRONGLY INTERACTING ELECTRONS 

In this paper we will be studying a strongly interacting 
one dimensional electron system coupled to either a ferro- 
magnet or a superconductor. Strong interactions imply 9 
that E a <C E p and therefore it is possible to be in a 
regime of temperatures that may not be small compared 
to the spin energy E a , but are still small compared to 
the charge energy E p . This means that we are no longer 
free to use standard bosonization procedures'^ for the 
electron operator because these formulas rely on the as- 
sumption that the relevant energy scales on which the 
system is probed are small compared to both spin and 
charge energies. 

For a strongly interacting system we require a formal- 
ism that applies to all energy scales relative to the spin 
energy (and in particular to fcgT > E a ), but only low 
energy scales to relative to the charge energy. Such a 
formalism was recently developed by Matveev, Furusaki, 
and Glazmanfi^ and applied to the evaluation of the spec- 
tral function and tunneling into a strongly interacting 
electron system at arbitrary temperatures relative to the 
spin energy^i In this section we review the essential ele- 
ments of their work and show how to represent the elec- 
tron operator in a strongly interacting one dimensional 
electron system. The general expressions for electron cor- 
relation functions obtained in this way provide a launch- 
ing point for numerical studies in the intermediate (with 
respect to spin) temperature regime fc^T E a , which is 
not well handled by existing analytical methods. 

A. Basic assumptions and general expressions 

As is typical of interacting one dimensional electron 
systems, we assume the Hamiltonian in the strongly in- 
teracting case is spin-charge separated;*^! H = H p + H a . 
Here, H p is identical to that of the LL, and is given by 
Eq. (TT|). On the other hand, the spin Hamiltonian at 
arbitrary temperatures is to a very good approximation 
given by a nearest neighbor antifcrromagnetic Heisenberg 
spin chain ) 14 ' 18 i 19 

H a = ^ J Si ■ Si + i, (2) 

where evidently the spin energy is set by J: E a = J. The 
basic idea is to represent the electron operator as a prod- 
uct of operators that describe the holons ^>(x) (spinlcss 
fermions that naturally arise in the context of strongly 
interacting fermions and the spin-incoherent regim o 20 ! 21 ) 
and the spin degrees of freedom Si . The holon operators 
(denoted by & , 4") by construction satisfy the equation 

&(x)*(x) = v4(aOVr(aO + i>\(x)i>i(x), (3) 

where ip s is the electron annihilation operator for elec- 
trons of spin projection s, and "01 the corresponding elec- 
tron creation operator. 
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The issue of how to bosonize the electron operator for a 
strongly interacting system earlier arose in the context of 
the large U limit of the one dimensional Hubbard model. 
Penc et al£2- wrote the electron creation operator as 



^1(0) = 4 s *t( ), 



(4) 



where Zq creates a site on the spin chain @ with spin 
projection s. The expression Q can be physically mo- 
tivated as follows. From ([3]) it is clear that the creation 
of an electron is also accompanied by the creation of a 
holon. However, electrons also carry spin so there must 
be a component of the electron operator that also cre- 
ates spin. This is accomplished by Z\ s . In general, one 

has Z\ s as the object that adds a new site to the spin 
chain between I — 1 and I. While this appears physically 
intuitive, the expression suffers from the drawback that 
it does not naturally account for the variation of electron 
density with position in a real electron gas^ Matveev, 
Furusaki, and Glazman showed^ that a remedy for this 
issue is to define the position at which the spin site is 
added to the chain @ in terms of the number of holons 
to the left of the site, 



l(x) 



tf(y)*(y)dy. 



(5) 



In terms of (J5j) the electron creation and annihilation 
operators are defined as 



ip s (x) = ^{x)Z u?0 ) < 



(6) 
(7) 



The operators given above explicitly account for the fact 
that the spins are attached to electrons, and the formu- 
las are valid at all energy scales. It is perhaps worth 
noting that even though the Hamiltonian is spin-charge 
separated, the electron operators are not written as a 
product of a spin piece and a charge piece because the 
"spin" pieces Z^ x ^ s also depend on the electron density 
via ©. 



B. Bosonizing the holon operators 

In writing Eqs.([H|),([7]) no assumptions have been made 
about the energy scale relative to the spin and charge 
energies. We now restrict our considerations to energies 
small compared to E p , but arbitrary with respect to E a . 
In this case, we are free to bosonize the holon sector. The 
Hamiltonian for the holons must then necessarily take the 
low energy forrr>2£ 



Hn 



dx 
2^ 



(8) 



where the interaction parameter K of the holon is related 
to the interaction parameter of the charge sector as 2 ^ 



K = 2K p , and the spinless fields 9 and 4> can be related 
to the holon density a o 21 i 24 



¥{x)*{x) = ~[k h F + d x e(x)\, 

7T 



(9) 



where the holon Fermi wave vector is twice the electron 
Fermi wavevector— k F = 2k F . The bosonic fields satisfy 
the commutation relations [0(x), d y <fi(y)] — iit5(x — y). 

Since we are interested in low energies with respect to 
the charge energy, the electron operator may be expanded 
about the two holon Fermi points at ±k F , 



= ^ R (x) + *l(x), 



(10) 



where "J r(x) destroys an holon near the right Fermi point 
and l(x) destroys an electron near the left Fermi point. 
The left and right holon operaters are bosonized as 



\J2not c 



_ e -i<t>(x) e ±i[k%,x+9(x)] 



where a c is a short distance cut-off of order the inter- 
particle spacing a. Combining the results of (J5|), (O, 
© , (flU)) and (jlip one obtains the bosonized form of the 
electron annihilation operator for spin s 



ip s (x) 



e -i<p(x) 



+6(x)} 



xZj 



l.s 



-i[kpX+9(x)] S 



l=^[k h F x+e{x)} 



(12) 



and an analogous expression for the electron creation 
operator i^l{x). Expression (fT2|) . however, it not quite 
complete as it does not account for the discreteness of 
the charge of the electron. This can be accomplished by 
interpreting 



Z, 



l = ±[k*x+6(x)} 



I 



Zi iS s -[kpx + e(x)] - 1 



(13) 

after which the full electron annihilation operator (in- 
clude both left and right moving parts) becomes 16 



where 



e -i<p(x) 

y / 2na c 



^- Zs ( q )e^+i^+ e ^\ (14) 



2tt 



Zs{q) = Z > 



i, s e 



-iql 



(15) 



i=-c 



The expression for the electron annihilation operator (fTij) 
is the key result of the this section, obtained earlier by 
Matveev, Furusaki, and Glazman, who also showed that 
in the limit of small energies compared to E a the expres- 
sion correctly recovers the standard LL formulas for the 
electron annihilation operator^ With ([T4]) correlation 
functions involving electron operators can be expressed 
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in terms of the correlation functions of the holon and 
spin sectors at arbitrary temperatures with respect to 
E a , but small energies compared to E p . See Appendix lAl 
and Appendix [B] for examples of the evaluation of the 
single particle Green's function near different types of 
boundaries using the bosonization formula (114j) and the 
boundary conditions to be discussed next. 

C. Open and Andreev boundary conditions 

In this work we will be primarily interested in two 
types of boundary conditions on the electron operators: 
(1) open or "hard wall" boundary conditions appropri- 
ate for tunnel junctions and (2) Andreev boundary condi- 
tions appropriate for adiabatically contacted (no electron 
backscattering) superconductors. Other, more general, 
"intermediate" boundary conditions are possible though 
less generic^ To help orient the reader, in each of the 
two cases above we will briefly summarize the result ap- 
propriate for the LL regime for the purpose of drawing 
contrast with the strongly interacting, finite temperature 
SILL regime. Our discussion of the boundary conditions 
extends the results of Ref . [lj| • 

1. Open boundary conditions 

Open boundary conditions, or "hard wall" boundary 
conditions result in electron waves that are perfectly re- 
flected at the boundary. This implies there is no charge 
current or spin current through the boundary. For con- 
creteness, let us assume that our boundary is located at 
x = with the interacting one dimensional system living 
on x > 0. Then an electron traveling to the left with 
spin s will be reflected to a right moving electron with 
the same spin s 

V>L, s (0)=e- l "Vi?,,(0), (16) 

with a phase shift rj that depends on details of the bound- 
ary scattering potential. 25 From the standard LL bo- 
zonization formulas (in our convention) 

1>h,lA*) = -^=e-^e H ^ 3!+e '^\ (17) 

we see that (fT6|) implies that S (Q) — const, and therefore 
also that P (O) = const and ^(O) = const. In computing 
correlation functions at the boundary, this means that we 
must take 9 P (0) and 9 a (0) to be non- fluctuating quanti- 
ties, allowing only (f> p (Q) and </>o-(0) to fluctuate. 

In the strongly interacting case, the reflection condi- 
tion (flu) implies that 9(0) = const, where we have used 
(fT2")l . Since the spin site I is related to the 9 field via © 
and ([5]) we have I = at x = 0. This implies any corre- 
lation function involving electron operators evaluated at 
the boundary will depend only on Zi—q >s or its Hermitian 
conjugate, and the fluctuations of the field (f> at x = 0. 



We will apply these boundary conditions in a calculation 
of the single particle Green's function in Appendix |A1 

For completeness, below wc give the expansions of the 
bosonic fields for a system of finite length L with open 
boundary conditions at x — and x — £: 24 i 26 i 27 

= *E^sin(^)(6 m -6tj + ^), 

rn — 1 

oo / ~ 

«*) = E^-s(=)( 6ro + & L) + $! (18) 

m=l V 1 

where 9°(x) = fW, [b n ,b\ n ] = S nm , and [$,N] = i. 
For a semi-infinite system, we take L — > oo and the 
discrete sums over m become integrals over momentum 
q m = rmr/L. 

2. Andreev boundary conditions 

In a certain sense the Andreev limit is the opposite 
limit of "hard wall" scattering from a SC-M, or SC-LL, or 
SC-SILL interface. Whereas open boundary conditions 
imply an incident electron is perfectly reflected, Andreev 
boundary conditions imply that an electron is perfectly 
absorbed by the SC (with a concomitant reflected hole 
of the opposite spin)^ However, compared to the open 
boundary conditions (|16| the Andreev boundary condi- 
tions are more subtle as they involve an energy scale, 
the superconducting gap A, that under many circum- 
stances cannot be taken to be infinitely large relative to 
the energy scales of interest (such as fcsT, or an applied 
voltage) in the interacting one-dimensional system. One 
of the most important consequences of a finite A is an 
energy-dependent reflection coefficient, which ultimately 
leads to the proximity effect in the normal material^ In 
the context of interacting one dimensional systems, An- 
dreev boundary conditions have been discussed by a num- 
ber of authors^^^L^ 2 ^^^ The conclusion of these 
works is that when a left moving electron with spin s is 
reflected as a right moving hole with the opposite spin 
there is an energy dependent phase shift e lq ^ (propor- 
tional to the momentum difference q with respect to the 
Fermi point )22 that multiplies a factor— e %x that encodes 
the phase \ 01 the superconducting order parameter (as- 
sumed non-zero for x < 0)r^ 

V>l, s (0) = (-1) /W ie < *e < «tyi,_.(0), (19) 

where £ oc 1/A is the superconducting coherence length. 
The function f(s) = for s =| and f(s) — 1 for 
s =1. The boundary conditions p9[) are valid only at 
energy scales much smaller than A, which we will assume 
throughout this work. 

Applying the boundary conditions (fT9|) to the LL 
case (with bosonized electron operator below (|16p ) gives 
9 a oc 9 S — 6—s = const and <j) p oc <p s + <fi_ s = const. 
Therefore, we find that much like the situation of perfect 
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reflection there is no spin-current through the interface, 
but there is is a net charge current. Moreover, analysis 
of singlet superconductivity and spin density wave corre- 
lation functions using (|19[) in the LL regime shows that 
there are suppressed spin fluctuations near (distances less 
than £) the interface.— It is perhaps worth noting that if 
the SC-LL interface has a very weak electron backscat- 
tering, the interactions in the LL tend to renormalize the 
interface scatteringi 25 i 32 i 34 

For strongly interacting electrons the bound- 
ary condition (JTHJ) implies (f)(0) — const and 

EiZi,sS(e(o)/w-i) = £iSj-/(0(Q)Ar-O =*► 

Zi, s = Zj at the interface, where we have again used 
(fl~2|) and also (fT3j) . We show in Appendix [B] that we 
recover our earlier results^ for SC-SILL junctions in the 
Andrecv limit using this formalism. 

Finally, for completeness and for unification of our 
notation, below we give the expansions of the bosonic 
field o 15 i 30 ' 36 for Andreev boundary conditions on a sys- 
tem of length L (see Fig. [3]) with superconducting order 
parameter phase difference \ — Xi ~ X2- 



Or) 
(a) 




(b m + bt) + e°, 



$(a:) 



+$(x), (20) 
[b^bln] = 6 nm , and 



where 

[8°, J'} = i. The topological number J' = (Nf+N L )/2+l 
and the total spin of the system M = (Nj — iVj.)/2 must 
satisfy the constraint J' + M =even.— We have also in- 
cluded the proximity effects 33 via the length £ in the 
field expansions that appeared earlier in the general ex- 
pression (|19| for the boundary conditions appropriate to 
Andreev reflection.— 

Note that compared to the field expansions for open 
boundary conditions (|18[) . the expansions for Andreev 
boundary conditions (|20|) have 9 and tj> "switched" . This 
can be understood in simple physical terms: 9 and 4> 
are conjugate fields so if one is a constant, the other is 
strongly fluctuating. Thus, the "switching" of the fields 
comes from "opposite" nature of the two boundary condi- 
tions in the charge sector. For open boundary conditions 
there is no charge current through the boundary, while 
for Andreev boundary conditions the charge current on 
either side of the boundary is unchanged by its presence 
because there is no electron backscattering there. Had 
we been concerned with the bosonization of the spin sec- 
tor, we would have found that due to the absence of spin 
current through the boundary that sector would have ex- 
pansions appropriate for open boundary conditions sim- 
ilar to [p]) ,15 ,30,36 

Having spent the time to develop the formalism used 
in our calculations, we now turn our attention squarely to 
the physics of hybrid junctions involving spin-incoherent 



Luttinger liquids. We first discuss a junction consist- 
ing of a ferromagnet tunnel coupled to a spin-incoherent 
Luttinger liquid. 



III. FM-SILL TUNNEL JUNCTIONS 







FM 'vW 


SILL 







x=0 



x=L 



FIG. 1: Schematic of the model we study. A ferromagnet is 
tunnel coupled at x = to a spin-incoherent Luttinger liquid 
of length L. As the magnetization vector m precesses about 
an effective magnetic field H spin is pumped into the SILL. 
Spin accumulation in the SILL will lead to renormalization of 
the magnetization dynamics. 



A. Spin pumping 

We consider the set-up shown schematically in FigQ] 
whereby a ferromagnet is coupled via a tunnel juntion 
to a spin-incoherent Luttinger liquid. A similar situation 
was considered by Bena and Balents in the context of a 
FM-LL junction^ They found that when the magnetiza- 
tion m of the FM acquires a time dependence (perhaps by 
the application of external fields), spin current is pumped 
into the adjoined LL at the rate 

,r>. , _ dm , dm . . 

(!(*)> =-A x mx — +A 2 — , (21) 

where A\ and Ai are temperature dependent parameters 
depending on the interface tunneling and interactions in 
the LL. Below we show that an identical expression is 
obtained for a SILL and we derive explicit expressions for 
A\ and A2 in this case. We schematically indicate how 
the temperature dependence crosses over from the LL to 
the SILL case. We note that (f2"Tj) was earlier derived 
for a non-interacting metal attached to a FM. In this 
case, A\ and Ai are temperature independent to lowest 
orde n 39 i 40 i 41 Thus, the temperature dependence of A\ 
and Ai encodes information about the interactions and 
may also indicate whether the attached 1-d system is in 
the LL or SILL regime. We now turn to a derivation of 
these results. 

We assume that the ferromagnet is itinerant, such as 
Fe or Co, and can be described by a Stoner-type model 
in which there is a different density of states for spin-up 
and spin-down electrons. The different density of states 
can be absorbed into distinct effective tunneling matrix 
elements, t s , for spin- up and spin-down electrons 
This also implies that the local action 4 ^ describing the 
spin up and spin down electrons are identical to that of 
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non-interacting electrons 



'FM 



R X] 



--Id 



2tt 



ICK)| 5 



(22) 



where /? is the inverse temperature and ^(w„) are the 
bosonic fields describing local fluctuations at the tunnel- 
ing point, x — 0, of spin up and spin down electrons in 
the FM. For the SILL we have the following spin-charge 
separated form of the local action 



SILL' 



(23) 



where < K p < 1 is the Luttinger parameter of the 
charge sector of the SILL, cj) p = (0 T + i )/v / 2 is the 
bosonic charge field, and S'sill is the action for the spin 
sector which will actually play no role in the evaluation 
of correlation functions in the spin-incoherent regime as 
it effectively drops out leading to "super universal spin 
physics" .— 

For an arbitrary spin quantization axis relative to the 
magnetization fa it is useful to define 



u± = (l±m- ct)/2, 



(24) 



which projects the spin quantization axis onto the mag- 
netization direction fa. Here u± is a 2 x 2 matrix and a 
are the Pauli spin matricies. The tunneling Hamiltonian 
then takes the form 43 - 45 



Hi-, 



(25) 



where F (F^) annihilates/creates an electron in the FM 
and tjj WO annihilates/creates an electron in the SILL. 
The tunneling is assumed to occur at x = as shown in 
FigfJ] The 2x2 tunneling matrix W is then 



The tunneling Hamiltonian 
explicitly as 



(26) 



can be expressed more 



x Fli> 8 +m{t)-Y,^Ft^f^ sl +h.c., (27) 



where fa{t) is the time-dependent magnetization of the 
ferromagnet, and u\ = (t + + t-)/2 and u 2 = (t + — 
t-)/\fh\. (From here onwards we will assume \m\ = 1.) 
The spin current operator is obtained from the relation 

*W = # = ff tun], where M = \^ s cr ss ^ s , is the 

spin density in the SILL at the boundary and summation 
over repeated spin indices is understood. It follows then 
that the spin current operator is given by 



I(t) 



— F± cr ss > ip s > H — —mF< ip B +—F<mx a ss , ip s > +h.c. 



(28) 



The spin current, = -ifdt'0(t - 

*0 ([/(*)> -^tun (*')]) j i s obtained 3 ^ from second-order 
time dependent perturbation theory which gives 

(T(t)) = m(t)Im(u* Ul )Re{C(0)} -J yC(w) e - ui x 

Im(u^ui)m(w) + -^-m(w) x fn(t)) (29) 



where C(lo) is the Fourier transform of the re- 
tarded Greens function C{t — t') = —iQ(t — 
t'){[F%{t)ilj s {t),il)l(t')F s (t')]). If we assume the typical 
frequencies of the magnetization precession are in the 
GHz range 4 - (which corresponds to energies of roughly 
100 mK) or smaller, then this is a small energy scale in 
the problem and we may safely expand C(u>) for small 
lu. It can be easily checked that the zero frequency terms 
cancel exactly, leaving only the contributions linear in lu, 
provided we drop terms proportional to lu 2 and all higher 
powers. Upon integration over frequency, the linear fre- 
quency terms are converted to time derivatives yielding 
expression (|2"Tj) where the coefficients A\ and A 2 are given 
by 



Ax 



A 2 = 



-C"(0)lm(^tti). (30) 



Since only the imaginary part of C(lu) is odd with respect 
to frequency, only this piece will contribute to A%, A 2 and 
those parameters will therefore be real quantities. The 
final step is to compute the temperature dependence of 
C (lu = 0) in the spin- incoherent regime. 

Since we are working within linear response, the com- 
mutator in C{t) = -iQ(t)([F^(t)i) s (t),ipl(Q)F s (0)]) 
is evaluated in the state where there is no 
tunneling between the ferromagnet and the 
SILL, and a hard- wall at i = 0. By writ- 
ing C(t) -ze(i)[(Ft(i)F s (0))(^(i)^(0)) - 

(F s (0)FJ(t))(ipl(0)ip s (t))], and noting that the long 
time behavior of both terms has the same functional 
dependence, we can easily extract the long time 
behavior— of C(t) using (j2"2")l and (|2"3"|) . (A more care- 
ful calculation that yields the same result is given 
in Appendix [X]) At finite temperatures (small with 
respect to the Fermi energy of the ferromagnet and 
charge energy of the SILL, but large compared to 

(* kBT 'V an d 



(nk B T/H) 



sinh(7rk B Tt/h) 
1/(2K P ) 

The crucial 



Hvp/L), we have 

(V'.(t)V'J(o)) ~ [ sinh ^ kB Tt/h) 

difference with the Luttinger liquid is that the exponent 
for the correlations (ip s (t)ipl (0)) have changed: in the 
case of a LL 1/{2K P ) is replaced by (1/K p + l)/2. The 
remainder of the computations carry through exactly as 
they would for a LL and we find for fku <C ksT 

Im[CV)] oc huj (k B T) 5SILL , (31) 
where S SILL = jr— 1 - This implies C"(0) cx T s ^ quite 
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generally so that 

Al{xT 5(T)M ; A 2C xT 5 ^lm(u* 2 u 1 ), (32) 

where S(T) interpolates between the LL and SILL 
regimes as the temperature is swept through the spin 
energy E a w J. The temperature dependence of the ex- 
ponent S(T) is shown schematically in Fig. [2] 



spin-coherent 

spin-incoherent 




5(T) 

(l/KpHl)/2-l 
1/(219-1 



FIG. 2: Temperature dependence of the exponent 8 appear- 
ing in the parameters (|32[) that describe the pumped spin 
current ()2ip due to time-dependent motion of a magnetiza- 
tion vector near an interacting acting one-dimensional system 
at finite temperature. For temperatures T less than spin ex- 
change J = E a , the characteristic spin pumping temperature 
dependence in the spin-incoherent regime crosses over to the 
Luttinger liquid form computed in Ref. [381] . Note that in the 
SILL, S(T) < for K p > 1/2 and the qualitative temper- 
ature dependence of the pumped spin current is remarkably 
different from a LL. 

It is interesting to note that 5 > for < K p < 1 
in the Luttinger liquid regime implying that less spin 
current is pumped as the temperature decreases. On 
the other hand, so long as the system remains in the 
spin-incoherent regime the opposite behavior may be ob- 
tained if K p > 1/2: since S SILL < more spin current 
is pumped as the temperature is lowered (for T > J). 
This is related to the diverging density of states at the 
boundary ^ 24 ! 46 when K p > 1/2. In gated cleaved-edge 
overgrowth quantum wires it appears possible to lower 
K p down to values of order 1/3 and so it may be pos- 
sible in the spin-incoherent regime to tune between the 
positive and negative exponent regimes £ This qualitative 
difference should be easily seen in experiment. 



B. Renormalization of Magnetization Dynamics 

Having computed the spin current pumped into the 
SILL by a time-dependent magnetization vector, it is im- 
portant to ask how the spin accumulation in the SILL 
in turn affects the magnetization motion. We address 
this question by computing the renormalization of the 
Gilbert damping constant due to the spin flow into the 
SILL. We again closely follow the notation of Bena and 
Balents 3 ^ to clearly establish a connection to the LL case. 
The Landau-Lifshitz- Gilbert equation for magnetization 
to precessing around effective magnetic field H is 



dm 
~~dt 



-7m x H + a^fh x 



dm 



where the spin current I = Iq — If, flows (Iq is given by 
(l2"Tj) for If, = and is non-zero only for time dependent 
to) into the SILL and M s is the saturation magnetiza- 
tion of the ferromagnet. Here 7 is the gyromagnetic ra- 
tio which is typically equal to its free electron valued 
7 = 2fj,B/h, in transition metal ferromagnets, and ao is 
the dimensionless Gilbert damping parameter in the ab- 
sence of the spin current. Its value is typically of order 
10~ 2 . The spin backnow due to spin accumulation in 
the interacting 1-d system can be described by boundary 
conditions on left and right moving spin currents and is 
expressed a s 43 i 45 



Ms 



iS(T) 



47T 



-fl s X TO, 



(34) 



where fl s is the spin chemical potential in the wire related 
to the magnetization by fi s (%) = M(x)/x(T), with x(T) 
the (generally temperature dependent) spin susceptibil- 
ity and ivexch describes the effective exchange coupling 
between electrons in the SILL and the FM4£ The current 
Ig(r)(n s ,T) arises from the electron tunneling contribu- 
tion to the spin current^ Note that the back-scattered 
spin current (134|) contains two terms: (i) a term arising 
from electron transfer from the FM to the interacting 
1-d system and (ii) a second term arising from exchange 
between the local spin density and the precessing magne- 
tization vector to. Since I = Iq — If, this implies the spin 
current itself has a contribution due purely to exchange 
effects, even in the absence of electron transfer^ 
Our goal in this section is to express (fg3"j) as 



dm . _ ~ _ dm 

= —7 m x H + am x — 

dt dt 



-7T> (35) 



7 



dt M s 



(33) 



and determine the renormalized parameters 7' and a. 
From (|3"3"|) it is evident that we must find I = I — lb- 
We have already computed Io, Eq. (|2"Tj) . in the previous 
subsection. Now we must determine h, ([34]) . which de- 
pends on /2 S and Is(t) (a*s, T). We start with fl s which 
is a function of the spin dynamics in the interacting 1-d 
system. 

It has earlier been shown that SU(2) invariant elec- 
tron backscattering leads to diffusive spin behavior in 
the weakly interacting regime i 43 i 49 In the strongly in- 
teracting regime where the spin sector may effectively 
be described by a Heisenberg spin chain coupled to 
phonon distortion o 9 ' 14 i 16 ' 50 the spin dynamics has also 
been shown to be diffusive^ In both limits, the mean 
free path I ~ v a /T. Following the argument o 20 ' 52 that 
the spin-incoherent regime can be approached from "be- 
low" (T < J), we expect the diffusion length I to saturate 
at the interparticle spacing, a, for T > J. Since the dif- 
fusion constant D s ~ lv a , this implies the diffusion con- 
stant in the SILL becomes independent of temperature 
and takes the value D s w av a . Note that in the limit 
v a — > 0, there is no spin diffusion, as the spin excita- 
tions cannot propagate in the system. Spin will simply 
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pile up at the boundary of the FM without moving fur- 
ther into the interacting 1-d system. Unlike our earlier 
results for the temperature dependence of pumped spin 
into the SILL, Eq.® with S SILL = 1/(2K P ) - 1, which 
were valid for vanishing spin velocity, here our results are 
qualitatively dependent on keeping v a finite, although 
still small enough to be in the spin-incoherent regime. 

We also assume the finite length SILL is characterized 
by a spin-flip time T s t (due to impurities, spin-orbit ef- 
fects, etc). The diffusion equation for spin in the SILL is 
then 



iup, s {x) = D s d 2 x fl s {x) - T s j-p, s {x), 



I at 



with the boundary conditions d x ft s = — y D X ( T ) 

x = and vanishing spin current d x fl s (x) — at x = L. 
See Fig. Q] for the set-up. The solution to this equation 
is simple to obtain and is given by 



/2 s (cc) 



where k = 



cosh[n(L — x) 



Ksinh[ K i] \D s x(T) 



1 



I(x = 0), (37) 



D at . If the precession frequency ui is 

small compared to the inverse spin relaxation time, then 
to a good approximation n w 1/ y/D s r s f. At the bound- 
ary x = we have fl s = £1 where £ = coth[KL] kD ^ t ^ ■ 

Note that d x jl s (x) ~ e~ KX T(x = 0) w e" 



T »fl(x = 

0) so the spin current decays exponentially with distance 
into the SILL on a length scale set by the product of the 
diffusion constant and the spin-relaxation time. For a 
long spin-relaxation time, this length scale can be large 
compared to the inter-particle spacing. 

We have now determined all parts of the spin back flow 
(|34|) . except the tunneling current I$(t)(h s ,T), which we 
now do. The tunneling current is proportional to the 
imaginary part of the Fourier transformed correlation 
functio n 43 ' 45 C(ui) defined below Eq. 



J, ( r)(A*.,T)«|ui| 2 Im[C(/i./2)] 
ocK|^ B T)Wsinh(^) 

S(T) 



r i 



oc n s \ Ul \ 2 (k B T) s W 



47rfc B T 



(38) 



where we have used the results of Appendix [X] and taken 
the limit [i s <C ksT. Evidently, the main effect of the 
spin-incoherent physics is to change the exponent S(T) to 
the spin-incoherent value, S SILL = l/(2K p ) — 1, so that 
we have 

I s(T) (jjl„T) oc f i s \u 1 \ 2 (k B T) 1 ^ 2K ^- 1 (39) 

in the spin-incoherent regime. Let us define T = 
Ig(T){lJ,siT)/ /j s tx (ksT) s which is temperature de- 



pendent and n s independent. At this point the deter- 
mination of 7' and a is identical to the LL case<^ We 
find 



7 



7 



and 



1 + 7 (Bi^2 - B 2 A 1 )/M S 



a a + 7 (BiAi + B 2 A 2 )/M S 



1 + 7(#iA 2 -B 2 A X )/M t 
where B t = (1 + £T)/[(1 + fT) a 



(40) 



(41) 



(^e X ch) 2 /(167r 2 )], 



(36) B 2 = (eX OX ch)/[(l+eT) 2 + (^c X ch) 2 /(16^ 2 )], and^,^ 
are given in Eqs. (|3QI32p . The temperature dependence in 
Bi, B 2 is entirely contained in T and £. See TableUbelow 
for a comparison of the LL and SILL regimes. 



TABLE I: Temperature dependence of key quantities in LL- 
FM and SILL-FM hybrid systems. Here c' is a constant that 
depends on L, r s /, and D s . 





<5 


D 3 


X 




T,Ai,A 2 


SILL 


-J- - 1 

2Kp 1 


const 


1 
T 


T 


ti'JLL 


LL 




1 

T 


const 


coth(c'v / T)v / T 


T s LL 



As with the case of a LL, for a SILL we expect there to 
be little renormalization of 7' relative to 7, but the small- 
ness of the Gilbert damping ao means this may obtain a 
significant temperature dependent correction depending 
on 5 SILL . 



C. Other FM hybrid structures involving a SILL 

One could easily imagine other scenarios such as a FM- 
SILL-FM junction, a FM-SILL-M junction, or even a 
FM-SILL-SC junction. However, because the spin trans- 
port is generally diffusive, for junctions whose length L 
is long compared to ^D s r s f , the two contacts to the 
SILL will more or less behave independently from the 
point of view of spin transport. To the extent that the 
two leads are coupled, it is evident from the discussion of 
the previous section that most of the physics (and corre- 
sponding general formulas) present for a LL system also 
apply to the SILL system only with some modifications in 
the temperature dependence of parameters appearing in 
the theory. While this may seem like a somewhat trivial 
result, it is not. The temperature dependence serves as 
a means to determine whether spin-incoherent physics is 
present in the system. In particular, we note quite gener- 
ally that if 1/2 < K p < 1 the temperature dependence of 
many quantities (e.g., T,A\,A 2 and those derived from 
them) change qualitatively relative to the expectations for 
a LL. Such qualitative differences should be observable 
in experiment. 

It is worthwhile to step back and emphasize some of 
the essential differences between the SILL and LL cases. 
First, we note that the different temperature dependence 
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can be traced to three quantities: (1) The correlation 
function C(t) defined below Eq. (gjj). (2) The temper- 
ature dependence of the spin susceptibility \- (3) The 
temperature dependence of the spin diffusion constant 
D s . So long as the SILL is tunnel contacted to the FM, 
the correlation function C (t) will always appear at lowest 
order in perturbation theory and carry along with it the 
characteristic temperature dependence of the SILL. This 
quantity appears in both the description of the pumped 
spin current and the magnetization dynamics. On the 
other hand, when we are specifically interested in how the 
spin propagates in the SILL [as we saw for the backscat- 
tered spin current ([34]) ] the diffusion constant D s enters, 
and also the spin susceptibility via the Einstein relation 
for the spin conductivity, a s = xD s - We remark that 
the spin transport is generically diffusive in the spin- 
incoherent regime, while it may be either ballistic or dif- 
fusive in the LL regime . 43 ' 49 

Having now flushed out what we feel are the central 
considerations and results for FM-SILL hybrid struc- 
tures, we now turn our attention to SC-SILL hybrid 
structures. 



feet Andreev reflection with weak (perturbative) electron 
backscattering at the SC-SILL interface. In the LL con- 
text, this was addressed by Vishveshwara et al~ in a 
SC-LL-M system who found dips in the conductance of 
the junction as a function of gate voltage. These dips 
correspond to multiple electron reflections off the SC in- 
terface and occur at voltages corresponding to integer 
values, eV — nhit/r, of the inverse electron traversal 
time, t = 2L/v, where n is an integer counting the num- 
ber of traversals in the LLi^I For a SC-SILL-M system 
we would expect similar qualitative effects, although with 
an overall suppression in the conductance. 13 ' 14 ' 15 Finally, 
we briefly remark that the conductance per channel of 
an adiabatic SC-SILL junction should drop to 1/2 the 
value expected for a SC-LL junction ; 33 ' 35 Gsc-Sill = 
^(2Kp)^%-. Recall that the conductance of the adiabatic 
SC-LL interface Gsc-ll = (2K P )^- generalizes the re- 
sult for the adiabatic SC-M interface 5 ^ Gsc-m = 2^ 
realized for K p = 1. 

B. SC-SILL-SC Tunnel Junctions 



IV. SC-SILL HYBRID STRUCTURES 

Throughout this work, we assume that all energy scales 
are small compared to the superconducting gap, A, un- 
less explicitly stated otherwise, such as in the limit of a 
short junction (defined below). We are primarily inter- 
ested in the case where the superconductor is contacted 
via a tunnel junction with a spin-incoherent Luttinger 
liquid. 



A. SC-SILL Junctions 

Let us begin our discussion of SC-SILL hybrid systems 
by considering the simplest case: where the FM in Fig. [I] 
is replaced by a SC. Earlier we studied the properties 
of such a junction in the Andreev limit.— We found a 
number of remarkable properties, including a completely 
universal (independent of the charge and spin Hamilto- 
nians) tunneling density of states in the spin-incoherent 
regime. 15 In the opposite limit of a tunnel contact be- 
tween the SC and SILL, no proximity effects of the SC are 
felt other than those perturbative in the tunneling. 53 (In 
contrast to the Andreev limit where the behavior of the 
pair correlations and tunneling density of states follows 
directly from the form of the field expansions imposed by 
the boundary conditions f 1 ^) As the proximity effects are 
already weak in the Andreev limit (they decay exponen- 
tially fast with distance, unlike the power law decay in 
the IAM-) we do not pursue the very weak proximity ef- 
fects in the spin-incoherent regime in the limit of a weak 
tunneling SC-SILL junction. 

On the other hand, it is worthwhile to briefly discuss 
the behavior expected in opposite limit of nearly per- 



(a) 



(b) 



Xi 




*2 




SILL ^\ 







@ 

FIG. 3: Schematic of two ways a superconductor can contact a 
spin-incoherent Luttinger liquid of length L. The phase of the 
superconducting order parameter at the left contact is Xi an d 
at the right contact x%- We assume there is a tunnel barrier 
between the SC and SILL. In (a) the SILL is "end contacted" 
and the zero modes (finite length) of the SILL generally come 
into play, while in (b) the SILL is "bulk contacted" and the 
zero modes play no role since the system length is effectively 
infinite. 

We now turn our attention to the final topic of this 
work, SC-SILL-SC junctions of the types shown in Fig. [3] 
Our primary focus will be on the behavior of the ac and 
dc Josephson effects in such a hybrid structure. In the 
context of SC-LL-SC structures there have already been 
a number of theoretical works^ 3 ^^^^^^^^ and 
some experimental results! 62 ' 63 ' 64 ' 65 ' 66 

Earlier we studied the dc Josephson effect in a junction 
like that shown in Fig. [3ja) only with Andreev boundary 
conditions at the interfaces of the SILL with the SCs*i£ 
We found that in spite of the fact that the pair correla- 
tions decay exponentially with distance from the bound- 
ary of the SC, the Josephson critical current scaled in- 
versely with the length of the SILL and was reduced by 
a factor of 2 relative to either the SC-LL-SC or SC-M- 
SC case (which have identical critical currents that also 
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scale inversely with the length)] 25 ' 30 We remarked that 
this "robustness" is essentially related to the fact that 
the superconductor phase difference only couples to the 
charge degrees of freedom in the SILL, which remain 
completely coherent. This effect can also be taken to 
support the notion^ that in some sense superconductors 
are spin-charge separated^ Yet another perspective on 
this result can be obtained by noting that the Andreev 
boundary conditions on a system of finite length can be 
mapped onto a system with periodic boundary conditions 
of twice the length. 30 In this way, the Josephson response 
becomes equivalent to the problem of persistent currents 
for spinless electrons in a ring threaded by magnetic flux, 
with the superconductor phase difference being simply re- 
lated to the flux. Therefore, the dc Josephson current for 
Andreev boundary conditions is essentially given by the 
physics of spinless electrons which is why spin-incoherent 
physics does not manifest in a dramatic way. The factor 
of 2 reduction in the critical current and Josephson pe- 
riod are a direct consequence of the mapping of the SILL 
onto spinless electrons. 

We now turn our attention to the case of open bound- 
ary conditions at the interfaces between the SCs and 
SILL. In this case phase information must be carried be- 
tween the two superconductors by Cooper pair tunneling 
between them. Because the pairs must hop from one SC 
onto the SILL, then onto the other SC, the Josephson 
current will occur at 4th order in this tunneling process. 
Moreover, because the Cooper pair tunneling process cre- 
ates a pair of electrons locally in the SILL near the bar- 
rier, many Fourier modes of the electron operator will 
come into play. This is in contrast to the case of An- 
dreev boundary conditions where only the zero mode en- 
ters in the evaluation of Josephson current . 15 ' 30 This cru- 
cial difference will turn out to dramatically suppress the 
Josephson response relative to what one would have for 
a LL when the SILL is tunnel coupled to SCs. In short, 
we find the critical current decays as an exponential of 
the junction length (compared to a power of the inverse 
junction length for a LL at zero temperature) with the 
characteristic decay length the inter-particle spacing. 



1. DC Josephson Effect 

The SC-SILL-SC tunnel junctions shown in Fig. [3] are 
modeled by the following Hamiltonian, H = Hsi + Hs2 + 
Hsill + Ht where -Hsi/2 are s-wave BCS Hamiltonians 
for the superconductors and Hsill is the sum of Eqs. 
([5]) and (f5]). The tunneling Hamiltonian Ht is given by 

H T = J2 T ^ 1 (x = Q)^ s (x = 0) 

S 

+ T 2 tpg 2 (x — L)ip s (x = L) + h.c. (42) 

We assume the SCs support a phase difference of x = 
Xi ~ Xi an d have the same superconducting gap A. 



The Josephson current is obtained from the relation 

J=-2ek B T°^, (43) 
ox 

where — e is the electronic charge and Z is the parti- 
tion function. We employ imaginary time perturbation 
theory, where terms dependent on \ appear in fourth 
order,— 

InZ = J rfr 1 dT 2 dT3rfT 4 r l 2 (T 2 *) 2 Ft 1 (r 1 - t 2 ) x 

11(0, L;t 1 ,t 2 ,t 3 ,t 4: )Fs2{t3 - t 4 ) + 
h.c. + similar terms, 

(44) 

where 11(0, L; n, t 2 , T3, T4) is a two-particle (Cooperon) 
propagator in the one-dimensional interacting system. 
The "similar terms" account for time ordered permuta- 
tions and spin projections (there are 23 of them, 4! all 
together). The propagation of Cooper pairs in the super- 
conductors is described by the time-ordered anomalous 
Green function F$, 

F Si (r-r') ee (T r Vf (r)Vf(r')) 

(45) 

where the expectation value is taken with respect to H$i 
and the electron operators in (T45|) are evaluted at x = 
for S\ and x = L for S* 2 . 

Of the 24 terms (plus their Hermitian conjugates) that 
appear in the partition function, Eq. (|44j) . physical con- 
siderations aid us in choosing the most relevant ones de- 
pending on the length of the SC separation distance L. 
There are two important length scales in the junctions we 
study that give rise to corresponding time scales: The su- 
perconducting coherence length £ and the junction length 
We define "long" to mean L ^> £ and "short" to 
mean L <§; £. We now turn to a discussion of the general 
expression for 11(0, L; n, t 2 , T3, r 4 ) in the limit of long and 
short junctions. 

For long junctions the tunneling into the SILL is "fast" 
and the propagation of the Cooper pairs is "slow" . The 
relevant two-body propagator is 

n lon «(0,L;T 1 ,T 2 ,T3,r 4 ) = 

(^(0,T 1 )^_ s (0,T 2 )^ s (i^3)V' s t (i,T4)), 

(46) 

where ip s is the bosonized electron annihilation operator 
given by Eq. (TTJ} and the averaging is taken with respect 
to Hsill- We have |n - r 2 | oc £, |r 3 - r 4 | oc £, and 
\n ~ T 4\ ~ |r 2 - r 3 | oc L. 

In the opposite limit of a short junction L <C £, the 
propagation through the SILL is "fast" compared to the 
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"slow" tunneling of the Cooper pairs, so the two-particle 
propagator approximately separates into a product of two 
single-particle Green's functions^ 

n shOTt (0,L;T 1 ,T 2 ,T3,r 4 ) 
« (^(0,ri)^(i,7a))(^_.(0 > 75,)^l,(L ) r 4 )) 

«G s (i,T 1 -r 2 )G_ s (L,r 3 -T4), (47) 

where we have |ri — T2I oc £, |r3 — r 4 | oc £, and |n — w 
|t 2 — r 4 | oc L (note the re-ordering of electron operators 
relative to ([46])). 

i?nc? contacted SILL when L>^: We first consider the 
case of an "end-contacted" SILL, shown in FigJ3^,. For 



this geometry the expansions (| 18|) of the holon field op- 
erators are the appropriate ones. Our assumption that 
all energy scales are small compared to A (specifically 
ksT <C A, and Hv p /L -C A) enables us to approximate 
the anomalous Green's function (|45[) with a delta func- 
tion in time, i 7 !s i ( T — T ') ~ 7riV;(0)e lXi <5(T — r') where 
Ni(0) is the density of states of the SC Si at the Fermi 
energy in its normal state. This approximation forces 
the times the electrons tunnel into or out of the SILL to 
coincide, T\ — t 2 and t 3 = r 4 . Without loss of generality, 
we let T3 = and define r = Ti . Using the expression for 
the electron operator, Eq.fTU), the propagator (|4*5|) reads 




The calculation of the Cooper pair propagator is carried 
out in the same manner as that of the Green's function 
presented in Appendix [A] We first evaluate the expec- 
tation value of the spin chain site creation/annihilation 
operators, 

mL) = {z h , s Zi 2 ,- s Zl- s Zt, s ), (49) 

in the spin-incoherent regime. To this end it is convenient 
to fully exploit the symmetries and boundary conditions 
of the problem. We first note that the open boundary 
conditions at x = 0, L force 9 =const at both these 
points. As before, we take this constant to be zero so 
that all 9 fields effectively drop out of Eq. (|4"5|) . The mo- 
mentum integrals in (|48|) can then be done trivially to 
give delta functions on the sites I4 which sets l\ = 1% = 
and I3 = U = kpL/ir. Therefore, we have 

5(0, L) = (^o.^o,- s 4. W .,- s ^la, s ) ' ( 5 °) 

which in the spin-incoherent regime can be evaluated 
quite simply by appealing to the physics of the spins. 
Here S(0, L) measures the amplitude for two spins of op- 
posite orientation introduced at x = L to arrive at x = 0. 
Deep in the SILL regime where all the spins are random- 
ized and spin exchange is highly suppressed, this is just 
the probability of finding neighboring sites of the spin 
chain housing oppositely aligned spins, 

H(0,L)=(±J =er^\ (51) 



(48) 

I 

from which it follows that 

n^ g (0,L;r,0) = e-^^g-u^o.rjgW.o)^ (52) 

The remaining correlations over the charge degrees of 
freedom can be computed using the expansions (|18p and 
the identity e A+B = e A e B e -[A,B\/2 valid when ^ fl j 

commutes with both A and B. The final result is 

n long f0 L-t OW^ 1 " 2 f 3 \^ 

(53) 

where 77 <C 1 is a short distance cut-off that ensures ultra 
violet convergence of the integrals over the charge fluc- 
tuations in the wire. Up to unimportant phase factors 
and other overall multiplicative constants, the Josephson 
current in the wire is then 

X jf*(l+£^)* (M » 

where the integral has been cut off at times r = L/v p , 
the time it takes for charge to propagate between the two 
ends of the SILL. Compared to the LL case, the most im- 
portant difference is that the critical current scales as an 

exponential of the length, J^,™ 8 1 critical oc e $~ ln2 (j^j 

2 _ 1 

rather than a power law (^r) Kp in the LL regime^ 
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A measurement of the length dependence of the critical 
current then serves as clear signature of spin-incoherent 
physics (or lack thereof), provided one has some knowl- 
edge of the density to infer k F = 2kp = n/a. We also 
note that the superconducting phase difference x appears 
as an argument to the sin function, rather than the saw- 
tooth form we found for Andreev contacts^ The sin form 
follows from the assumption of tunneling contactsj 25 ' 57 i 58 

Bulk contacted SILL when L 3> £: We now consider 
the case of a "bulk-contacted" SILL, shown in FigJSjD. 
For this case, we assume that the perturbations induced 
by the contacts are irrelevant, which is so for K p > 1/2. 20 
Otherwise, if the perturbations induced by the contacts 
are relevant this case reduces to an "end-contacted" SILL 
discussed above. 

Starting again with we note that translational 
symmetry implies that the spin correlations satisfy 



(55) 



S(0,£) = (Z luS Z l2< _ s Zl_ s Zl 



z, 



z. 



zl 



l\— l4,s — — S 0, — s 0,s 



which allows us to shift the sums: h — h ~ h, h — h — h 
in (|48[) . The summation over 1% leads to a delta function 
setting (72 = 93 and the summation over Z4 leads to a delta 
function setting qi — (74. Making the change of variables 
q = {qi + <?2)/2 and q = qi — q^ one can do the integration 
over q which sets l\ = li- The result is 

xe- l2 ( 1+ -) fc ^ L (e l2 {( 1+ -) [e(o ^ ) " 9(L ' o)1 " [0(o ' T) " 0(i ' o)1 J'). 

(56) 

The expression (|56l) is valid quite generally under the 
assumption that the contacts are irrelevant perturbations 
and do not lead to the "end contacted" result described 
earlier. If we now specialize to the spin- incoherent case 
we have 



E(0,L) = {Z, s Z,_ X-X.> = (I 



(57) 



Substituting (|57[) into (|56[) and integrating over momen- 
tum, 



V 27ra c 



IS, p 



-5[l + ^(L + e(L,0)-9(0,T))} 

2)1 7T 

Xe i2{[e(o,T)-e(i,o)]-[0(o,r)-^(i,o)]}\ 



(58) 



Finally recalling that for L 3> a we can replace the sum 



over I by an integral; 2 ^ we find 



n£S(o,ii;r,o) 



/ 1 \ 2 1 *|£ 

V 2na c 



1 

r 



■ In 2 



x / e -^(e(L,0)-e(0,r))]\n2 e i2{[e(0,r)-9(LM-WO,r)-<l>(L,0)]}\ 



(59) 



which can be evaluated (up to phase factors) to give 



n 



long 
bulk 



(0,L;r,0) oc 



/ 1 



\ 2na c 



■ In 2 



L 2 + V 2 T 2 



iKr, 



where lKp ^^-^f 
rent in the wire is then 

Jl°:iix) cc GiG 2 e 

dx 



(60) 

The Josephson cur- 



2-/, 



sin(x) 



(61) 

where we have again cut off the integral at times 
t = L/vp. Evidently, the critical current again 

1 lk I critical 



scales as an exponential of the length, i7b°"] S| 



\ 2 1K„ 



rather than a power la\ 
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\l) Kp m * ne re gi me - This Luttinger liquid result 
corrects the result originally obtained by Fazio et aL 36 ^ 5G 
by taking into account the proximity effect. Comparing 
the critical currents of the end contacted (IS~H) case with 



the bulk contacted case (IBTj) , one sees that the length 
dependence is not materially different: they only differ 
in the inconsequential power law that multiplies the ex- 
ponential. 

End contacted SILL when L <C £: To compute the 
Josephson current we must evaluate the single-particle 
Green's functions that appear in (|Tf|). Since £ is a prop- 
erty of the SC which we assume is tunnel contacted to 
the SILL, it is independent of the properties of the SILL. 
We will further assume a <C L <C £. The opposite limit of 
the wire length being shorter than the inter-particle spac- 
ing is not well motivated physically. The single-particle 
Green's functions appearing in (|47|1 can readily be eval- 
uated following the method of Appendix [X] This gives 



^fend 



1 



2ira c 
1 



h L 



■ In 2 



a -i[0(r)-^(L)]\ 



>1 



2-KOLc" \\ + e -r) e -v P TTT/L 

which immediately leads to the Josephson current 



(62) 



dx 



■ In 2 



— ) sin(x) 

Otr. 



1 + e- r >e- x 



(63) 
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with a critical current J e tj rt I critical oc e - 2 ^ ln2 (j^f . 

Bulk contacted SILL when L <C £: For a bulk contacted 
SC in the short wire limit we again apply the formula 
(|4T|) where the Green's functions are those appropriate 
for an infinite system, computed earlier in the literature 
in Ref. @[2l|. The result is 

JhuVt (x) oc AGiG 2 e - j sm(x) 



tlx 



1 + x 2 



(64) 



which evidently has 



J -short I 
bulk I critical 



r 



1 

4A'„ 



K n 



the 

' ln2 (if) 



critical current 

2T K -1 

p , where 



1 



Comparing the results for bulk and end contacted, as 
well as short verses long wires, we find the most impor- 
tant result is how the Josephson critical current scales 
with the length of the wire. Suppressing the relatively 
unimportant power law decay that multiplies the domi- 
nant exponential decay we find 



jj-short 



critical OC & 



critical CX e 



In 2 



(65) 
(66) 



so that the decay of critical current with the length of the 
SILL is twice as fast as for a short wire. This originates 
in the fact that for a short wire electrons propagate in- 
dependently so the spin incoherence affects each electron 
independently, rather than a single coherent pair as in 
the case of the long wire. It is also worth emphasizing 
that the assumption of tunneling contacts always results 
in a Josephson current that is proportional to the prod- 
uct of the bare conductances of the two contacts with 
a sinusoidal dependence on the phase difference of the 
superconducting order parameter, J(x) oc GiG2sin(x). 
Therefore it is the length dependence of the critical cur- 
rent (|66p that provides a smoking gun signature of spin- 
incoherence in the dc Josephson effect with tunnel con- 
tacts. For Andreev (adiabatic) contacts, the length de- 
pendence is identical to those of a LL or non-interacting 
one dimensional system. 1 — Instead, the SILL physics is 
manifest clearly in the flux dependence of the Josephson 
current which takes on a saw-tooth form of half the usual 
period^ 



2. AC Josephson Effect 

The ac Josephson effect occurs when there is a finite 
voltage V across the SC-SILL-SC system. The Josephson 
phase acquires a time dependence x = 2e^ leading to 
a Josephson current that oscillates in time. A sub-gap 
dissipative current is also induced, but at small voltages 
this can be estimated to be small^ As in the case of 



the dc Josephson effect, the qualitative features of the ac 
Josephson effect also depend on whether the wire is short 
or long, as defined in the dc case. 

Long wire case: We have seen that the qualitative fea- 
tures of the dc Josephson current in the spin-incoherent 
regime for bulk and end contacts are not very differ- 
ent. That is also true of the ac Josephson effect. The 
ac Josephson current can be computed from 36 



J(t) = 4ir 2 ev] 



G1G2 

(4^)2 



xRc 



/>oo 

Y j ±e ±l2eVt dt'e* ieVt 'lL(t') 



(67) 



where IT(i) is the Cooper pair propagator, Eq. (|46|) . eval- 
uated at real times. The ac Josephson current can be de- 
composed into sinusoidal and cosinusoidal components, 



J(t) = ATr 2 ev] 



G1G2 

(4^)2 



[J s sin 2eVt + J c cos 2eVt] , (68) 



where 
J s = — Im 

J c = Re 



J dt' {e- ieVt 'u{t') + e +ieVt 'u(t') 
' (e- ieVt 'll(t') - e +leyt 'n(i')) 



dt 



*69) 
(70) 



In the Luttinger liquid regime J s and J c oscillate with 
v oltage ac ross the junction, as does the amplitude J a = 
\/J 2 + J 2 - 3 — The frequency of the oscillations of J a with 
voltage depends on the relative spin and charge veloci- 
ties. One period occurs when the spin and charge parts 
differ by 2?r, that is when eV = 2n/(L/v a - L/v p ). 36 
In the spin-incoherent regime, we have v a — > 0, imply- 
ing vanishing voltages will lead to oscillations and they 
may cease to be observed. The lack of amplitude oscilla- 
tions will persist into the spin-incoherent regime, which 
has effectively only one velocity, the charge velocity. On 
the other hand, if v a is not too different from v p (say, 
v a = Vp/10), then one can expect to find a tempera- 
ture dependence of the voltage oscillations that reveals 
spin-incoherent physics in a way analogous to Coulomb 
drag££ or charge fluctuation noisei^ - At temperatures be- 
low the spin energy, there will be amplitude oscillations 
in the Josephson current as a function of voltage, while 
for temperatures above the spin energy there will be no 
such oscillations because the spin mode effectively does 
not propagate. This "washing out" of the high frequency 
(because of the ratio v a jv p <C 1) oscillations with tem- 
perature is the signature of spin-incoherent physics in the 
ac Josephson effect. 

Short wire case: For a short wire, the voltage de- 
pendence is independent of the properties of the 1-d 
system^ be it a LL or a SILL. In this case the ac Joseph- 
son current is 



J(t) = -K(eV/2A)JJ0) sin(2eW), (71) 

7T 
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where J C {Q) is the zero voltage critical current for the 
short wire given in Sec. IIVB li and K(x) is an elliptical 

integral. Recall that for a short wire, J C (Q) oc e~ 2 ~^~ n2 . 
Thus, the ac Josephson effect is only effective at revealing 
spin-charge separation in the long wire limit, and there 
are no new spin-incoherent features that appear relative 
to the LL aside from the length dependence of the dc 
critical current that enters in (|7ip . 

V. DISCUSSION 

In this paper we have touched on what we believe are 
some of the most easily observed consequences of spin- 
incoherent behavior in ferromagnet/spin-incoherent Lut- 
tinger liquid and superconductor/spin-incoherent Lut- 
tinger liquid junctions. For the case of FM-SILL junc- 
tions, we computed the spin current pumped into the 
SILL as a result of magnetization dynamics and the ef- 
fect of spin accumulation in the SILL on the parameters 
governing the magnetization dynamics. We found that 
for interactions in the SILL with 1/2 < K p < 1 the 
temperature dependence of the spin current and magne- 
tization dynamics is qualitatively different from the case 
of a LL and should thus be observable experimentally. 
Some of the key differences between FM-SILL and FM- 
LL systems are summarized in Table HI The crossover 
from FM-LL to FM-SILL in the exponent S governing 
the temperature dependence of several key quantities is 
shown in Fig. [51 

In the case of SC-SILL junctions our results greatly ex- 
tend those obtained earlier by usJ£ In that earlier work 
we were concerned only with the case of adiabatic (An- 
dreev) contact of the spin-incoherent Luttinger liquid to 
the superconductor. Here we have developed those re- 
sults further and also discussed the opposite limit of tun- 
nel contacts to the SC. In the tunneling limit we have 
computed the ac and dc Josephson response in the ge- 
ometries shown in Fig [3] We find that in contrast to 
the case of adiabatic contacts, the tunnel contacts lead 
to a Josephson critical current that is exponentially sup- 
pressed with the length of the SILL region. This dif- 
ference arises from the fact that in the Andreev limit 
the dc Josephson effect is determined solely by the zero 
modes of the Hamiltonian, while for tunnel contacts the 
non-zero modes generally dominate the response. These 
non-zero modes enter because in the tunneling process 
an electron is created locally near the tunnel barrier and 
thus requires many wavevectors to build its wavepacket. 

With the aim of providing a general discussion of 
junctions of ferromagnets or superconductors with a 
strongly interacting one dimensional system, we have 
couched many of our calculations in the recently devel- 
oped scheme for bosonizing strongly correlated electron 
systems.* 6 This formalism is valid for arbitrary tempera- 
tures with respect to the spin energy E a , but requires 
that the temperatures remain small compared to the 
charge energy E p . In order to adapt that formalism 



to the systems we discussed here, we extended those re- 
sults to express open and Andreev boundary conditions 
in a language valid for arbitrary temperatures with re- 
spect to the spin sector. Through several examples, we 
showed how various correlations functions could be eval- 
uated and verified that in the spin-incoherent regime the 
results properly reduce to the results obtained using the 
world-line pictured 

One of our primary motivations for using the language 
of Ref. [l6[ is to help provide a clear starting point for 
numerical studies of the spin-incoherent regime and the 
many interesting (and probably experimentally relevant) 
crossovers that occur between it and the Luttinger liq- 
uid regime. Within this formalism the charge physics 
can be computed analytically via a standard bosoniza- 
tion scheme, but the spin sector must be addressed nu- 
merically for temperatures of order the spin energy E a . 
The types of correlation functions that must be com- 
puted numerically are those that add and remove a site 
(or multiple sites) from a spin chain, such as (Z/ l S Zj 2 s ), 
that appear in the evaluation of a single-particle (or 
multi-particle) Greens function. (See Eq. (|A4[) . for exam- 
ple.) We would like to emphasize that numerical studies 
of strongly interacting one-dimensional electron systems 
with appreciable temperature compared to the spin en- 
ergy is an entirely untouched area and is now ripe for 
investigation. 

On the experimental side, the lack of a clear experi- 
mental "smoking gun" observation of the SILL remains 
a key issue to be addressed. However, there are mount- 
ing experimental indications we are not far awa y 10 ' 11 ^ 2 
Our best numerical estimates*^ suggest that many low 
density quantum wires sit right on the edge of the spin- 
incoherent regime and perhaps all that is needed is a 
focused experimental effort to search for its signatures, 
rather than any key technical breakthough. One of our 
aims in this work it to highlight certain classes of hy- 
brid structures where spin-incoherent physics should be 
observable. 

Finally, we would like to close with what we regard as 
some of the outstanding theoretical issues surrounding 
the spin-incoherent Luttinger liquid. Perhaps the main 
one is the behavior on temperature scales ksT w E a that 
we already alluded to above. Related to this is a bet- 
ter understanding of the crossover between the Luttinger 
liquid and the spin-incoherent Luttinger liquid regimes. 
Both of these will likely require a numerical attack as 
there are no obvious analytical methods available to ad- 
dress them. There is also the issue of spin-orbit coupling 
that has so far received no attention. For very strong 
spin-orbit coupling is there novel behavior in the regime 
£„ « AbT < Eso,E p ? The subject of noise in hybrid 
structures involving a SILL will be discussed in a forth- 
coming workj^i 
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APPENDIX A: EVALUATION OF C(u) USING 
BOSONIZATION FOR STRONGLY 
INTERACTING ELECTRONS 

In this appendix we compute the correlator C(t) and 
its Fourier transform C(u>) using the general formalism 
for bosonization of strongly correlated electrons in one di- 
mension developed by Matveev, Furusaki, and Glazman 
who also applied it to the evaluation of the single parti- 
cle Green's function for an infinite system and its Fourier 
transform . 16 i 17 We summarized the main results of the 
bosonization scheme in Sec. [TTJ The calculation here is 
for a semi-infinite, or finite system with open boundary 
conditions, so is different in detail from what has been 
discussed in Refs. [313, 

but the basic elements of the 
bosonization are the same. This appendix is meant to 
illustrate clearly in a specific example which parts of cor- 
relation functions can be computed analytically at finite 
temperatures and which pieces in general require meth- 
ods yet to be developed, or a numerical attack. As was 
already discussed in Sec.[TT]the chief difficulty lies in com- 
puting the correlations at arbitrary temperatures in the 
spin sector. It is hoped that the details given here will 
provide a good starting point for those skilled in numerics 
to enter the study of strongly interacting one dimensional 
systems where there is currently no quantitative under- 
standing of the regime fc^T w E a -C E„. Given the 
typical values of T, E a , E p presenl&iiiiiL^ii 7 ^ in quantum 
wires, this "intermediate" temperature regime may turn 
out to be the most relevant experimentally. 

In the evaluation of the spin current (|21[) (pumped 
from a ferromagnet into an adjoined one dimensional sys- 
tem coupled via a tunnel junction), the correlation func- 
tion 

C(t) = -ie(t)^<[Ft(t)^(t),^(0)F s (0)]) (Al) 



arises in lowest (second) order perturbation theory. It 
contains the difference of the product of single particle 
Greens functions for both the FM and the interacting one 
dimensional system 

c( t ) = -ie(t)j2[{FKm(om s mt(o)) 

s 

-{F s (Q)FKt))(i>t(0)Mt))}, (A2) 



in Sec. ITTCTl The correlators (Fj(t)F s (0)) and 
(F 8 (Q)F$(t)) can be computed by standard bosonization 
methods. 1 The result is 



(Fl(t)F s (0))=i{F s (0)Fl(t)) 



i (irk B T/e c ) 
2ira c sinh(^f^) 



where e c = hvpja c is a high energy cut off. 

Our real objects of interest here are the boundary 
Greens functions G+(t) = (ip a (t)ipt(Q)) and Gj(t) = 
(-01 (O)"0s (*)} evaluated at arbitrary temperature with re- 
spect to E a , but small compared to E p . To evaluate these 
we make use of expression (|14p for the electron operator. 
Straight forward substitution gives 



^2 g-i(giii— q2h) 



x(e i[(1+ ^ )e(t) - 0(t)1 Z ;i s Zj e~* [(1+ ^ )9(0) ~< Wo)1 }.(A4) 



Open boundary conditions at x — prevent fluctuations 
of the 9 field which sets it to a constant, which we take 
to be 0. With 9 set to a constant, the integrals over 
momentum can be done trivially resulting in <5-functions 
for Zi,?2 which kill those sums and selects h = h = 0. 
The final result is remarkably simple 

Gf(t) = ^Z o At)zU0))(e- lWt) - m] )- (A5) 
Similar manipulations yield 



G-(t) 



(zlMZo,s(t))(e~ lWt) " m] ), (A6) 



where we have made explicit the time dependence in the 
spin correlators involving Zq s and Zo tS . The general the- 
oretical challenge is then to compute those correlation 
functions which involve adding and removing a site from 
the spin chain, which is an interacting many-body prob- 
lem analogous to Fermi edge physic s 73 i 74 i 75 only it in- 
volves knowledge of the Hamiltonian on potentially many 
energy scales. However, for E a <C E p the time evolution 
of the spin degrees of freedom are very slow compared to 
the charge degrees of freedom, regardless of temperature, 
and can be neglected for times < H/E a . Note that in the 
limit E(j — > 0, the spin dynamics can always be neglected. 
Furthermore, if fc^T ^> E a , then the spins become es- 
sentially non-interacting, and therefore non-dynamical. 
Then the correlation functions simplify considerably. In 
the spin-incoherent regime one ha a 17 i 21 ' 76 



and 



(Zo, s Z s ) - 1, 



(A7) 



(A8) 



where the brackets (•) denote a thermal average com- 
puted with the open boundary conditions described 



independent of time which can be straightforwardly 
generalized to include an applied external magnetic 
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fteld.^d&lZZL The result is 



(z ,.zL) = 1, 



(A9) 
(A10) 



where p s is the probability of having spin projection s. It 
takes the values y>f = 1 — p± = e -E Z /\ B T +1 i where Ez is 
the Zeeman energy. Note that the correlation functions 
(Zq 8 Zq ;S ) and (Zq iS Zq s ) in the spin-incoherent regime 
are identical to their values in the infinite system^ This 
is because in the spin-incoherent regime the spins are 
non-dynamical, so the open boundary conditions which 
tend to suppress fluctuations have essentially no effect on 
the spins which are rendered non-dynamical by fcgT 3> 
E a . 

Finally, we are left to evaluate (e - *^') - ^ )]), which is 
identical to (F^(t)F s (Q)) and {F„(Q)F%(t)) only with the 
exponent changed 



,-i[^(t)-0(O)] 



H) 



;\1/2K P 



2na c 



^k B T/e' c ) 
sinh(^) 



l/2K„ 



(AH) 

where e' c = hv p /a c is a high energy cut off for the charge 
sector in the SILL. 



Combining the results ([A!]) . (TAT)) . ([A"8]) . and (JXTTJ, we 
find 



C(t) 



i / 1 

-i6(t)(-i) 2ff p 



27ra f 



(irk B T/e c ) 
sinh(^p) 



J +2 



(A12) 



where (5 s/ii = l/2K p — 1 is given in Table Q] and 
e c is an effective high energy cut off given by e c — 
[ic(^' c ) 1 ^ 2Kp ] 2Kp ^ 2Kp+1 ^ _ p or making the following for- 
mulas more compact, we define 8 = 8 + 2. 

We are now ready to compute the Fourier transform 
of pOL2]) . C(w) = dte iuJt C(t). Making the change of 

variables X — F T t, we have 



This gives 



^ (2 ^-kUbt) 



r 1 







s_ - huj 

2 L 2irk B T 



I 5 Huj 
T \2 +l 2^f 



I 8 hw 
M 2 + *2^T 



(A15) 



where the sin can be expanded to pick out the real and 
imaginary parts. Selecting the imaginary part of C{lo) 
we find 



Im[CH] oc 



1 



(27ra c ) 2 \nk B T 



irk B T \' 



x sinh I — — — ) 
\2irk B Tj 



8 hw \ 
2 +l 2nk B T J 



which for too <C k B T gives 

hn[C(w)] oc hw(k B Tf 



(A16) 



(A17) 



in agreement with our (f3Tj) . and the general analytical 
form earlier obtained by Bena and Balents 3? up to the 
exponent describing the temperature dependence which 
is changed in the spin-incoherent regime, as illustrated in 
Fig! 

The calculation above is valid provided the tempera- 
ture is larger than the level spacing of the finite length 
wire: k B T > hv /L. For k B T < hv p /L the finite length 
results of Ref. [26j can readily be generalized to finite 
temperatures. No simple power laws emerge, but rather 
a more complicated dependence involving hyperbolic trig 
functions. 



APPENDIX B: SC-SILL JUNCTIONS IN THE 
ANDREEV LIMIT 

In this appendix we show that the bosonization scheme 
introduced in Sec. [TT] recovers our earlier results for the 
single-particle Greens function, the tunneling density of 
states, and the pair correlations^ 



CM 



(27ra c ) 2 



irk B T 



■nk B T 



dX- 



The integral is standard 



[sinhpf)^ 
(A13) 



dX- 



= 2 s - 1 T{\-8) 



[sinh(X)]* 



s_ ■ huj 

2 L 2-Kk B T 



r f 1 — - — i hu) 

1,1 2 l 2-nk B T 



The result of the integral can be transformed to a more 
convenient form using the following identities for complex 
number z: tt = sin(7rz)r(z)r(l — z) and T{z)* = T(z*). 



1. Single- particle Greens function 

We are interested in computing the Fourier trans- 
form of the single-particle Green's function Gf(x, x',t) = 
(tps{x, t)ij}\ {x 1 , 0)) where here x and x' measure the dis- 
tance from the boundary with the SC. Taking x = x' and 
using the general formula (I14|) we have 



G+(x,t) 



2ira r 



dqi 
2tt 



2ir 



-i(qih—q2te) 



x(Z h S Z\ )( e *[(i+^)[4z+0OM)]-0OM)] 

-i[(i+f-)[fcj.x+e(x,o)]-0(x,o)]\ (- B1 ^ 



xe 
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where the 9 and qb fields have the expansion given in ([20]) . 
We take the limit L — > oo as we did in Ref. [15[ so that 
the zero modes play no role. Since the holon sector is 
described the the Gaussian theory ([5]) we can make use 
of the identity (e lA ) = er^ A ^ 2 for operator A and evalu- 
ate the correlation functions at zero temperature with re- 
spect to the charge energy but infinite temperature with 
respect to the spin energy: E p ^> ksT — > ^> E a ^ 0, 
that is we take both the spin energy and the temperature 
to zero but alway maintain kgT 3> E a .— In the spin- 
incoherent regime, the spin correlations are translation- 



ally invariant so we have (Zi ls Zj^ s ) = {Z^_ 

l/2' ll ~ l2 '. Making a change of variables l = l\ — I2, the 
summation over the sites of the spin chain results in a 
delta function, 2ir5(qi —52), which immediately kills one 
of the momentum integrals and sets q\ = 02- We re-label 
the remaining momentum variable a, 



At the boundary, we showed in Sec. Ill C 21 that Andreev 
reflection implies that Z/.j = Zj Since I = 9(x = 0) at 
the boundary we may fluctuate, we expect also that the 
result Z/.j = zj 1 approximately holds near (within £) of 
the boundary. Hence, 



F(x) 



2na r 



dqi 
2ir 



dq2 --i(qih-q 2 l2) 



00 2n 11 



x (e 



— i<p{x) 



Z hA e 



i[(i+^)[k h F x+e(x)\ 



xe-^zl^ 1+ ^ kh F* +e ^}.(B5) 



Gf(x,t) 



2TTa c J_ 00 2n ^ 

L — — QC 



\l\ e -iql 



(B2) 



Next, we define a variable Y = 1 + q/ir and re-express 
the integration in terms of this variable, 



We now assume that we are in the spin-incoherent regime 
where the spin correlations are translationally invari- 
ant: {Z luS Z} 2 J = {Z h ^ s zl s ) = 1/2^-^1. Then 
repeating the steps used above to compute the single 
particle Greens function (changing the summation vari- 
able I = Zi — Z2, killing the momentum integral with 6- 
function), we find 



G+(x,t) = 



2na r 



I -ilnY 



xe 



f f i <*-«(-!, 

Z l=-oo 



i[Ye(x,t)-4>(x,t)] -i[Y8(x,0)-(p(x,0)]\ 



(B3) 



Finally, noting that 6(wl - [9(x,t) - 9(x,0)}) = 
^ e -*[*i-(0(x,t)-«(x,o))]y we see that Eq _ ,[b1J) cor- 
rectly reproduces Eq.(4) of Ref. [HI after we recall the 
relations (j) p /V2 = 4> and V29 p = 9. 9 - 20 The identical 
result therefore follows for the Greens function and the 
tunneling density of states derived from it. 



Pair correlations 



F(x) 



1 f°° ^1 V(e-^2-l'l 
2na c J_ 00 2n ^ 



xe 



-OO -■■ ; 

i2l[k^x+6(x)] -i2<j>{x)\ 



(B6) 



We then make use of the integral representation 



of the (5-function, f™^ |f Y,i 
6(1- ^[kpx + 9(x)]) to obtain 



,-iq(l-%[k%x+6(x)]) = 



Here we compute the pair correlation, F (x) = 
— (ip-\(x)?pi(x)), a distance x from the superconductor 
boundary. By direct substitution of (fT4")) and (fTS)) we 
have 



F(x) 



1 



2ira r 



dqi 
2ir 



dq : 



2tt ^ 

h,h 

)[k%x+6(x)]-<t>(x)] 



xe 



-i(fi(x) 



Z h,l e 



h,1 e 

i[{l+^-)[k%x+$(x)]-4>(x)] 



.(B4) 



F{X) W ~ 2^~ E< 2 ~ ( 1 " zV& x + « 

c 1 



which for x S> a, a c can be evaluated by taking the dis- 
crete sum to an integr al: Y\ — > J dl. Doing so we recover 
the result (9) of Ref. [l5| after again recalling the rela- 
tions (j) p /y/2 = <f> and V29 p = 6>. 9 - 20 
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